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Research  Area 


Mode  Competition  and  Multi-Frequency  effects  in 
high  power  fast  wave  devices 


Specific  Applications: 

1.  High  power  gyrotrons  for  plasma  heating 

2.  Parametric  Mode  interactions  in  harmonic  gyrotrons 

3.  Multi-frequency  effects  in  fel  oscillators 

Modeling  of  harmonic  gyroklystrons,  phase  locked  gyro 
oscillators,  and  gyro-twystrons 


4. 


NONLINEAR  SIMULATION  OF 
HIGH  POWER  GYROTRONS 


RESEARCH  TOOL: 


‘Time  dependent  multi-mode 


V  W  Vi  V 


CODE  FEATURES: 

*Axiai  RF  field  profile  determined  seif-consistentlv 
*Time  dependent 

^Multi-Mode  (frequencies  are  equally  spaced) 

^Realistic  boundary  conditions 
*AC  space  charge 

*DC  space  charge  (voltage  depression) 

^Energy  and  Pitch  angle  distributions  for  incoming  beam 
*Window  reflection  and  delay 
*Non  zero  beam  rise  time 

*Mode  coupling  due  to  non  uniform  wail  radius 


THEORETICAL  APPROACH 


•  Radiation  field  is  expanded  in  transverse  modes  TE,-., 

C  P 


ip  &  ip 


Epp(z.t)exp(-iLuept)ezxV(Jc(kipr)e 


•  Particle  trajectories  are  followed  in  prescribed  fields, 
AC  beam  current  calculated 


•  Axially  dependent  field  profile  E^z^t)  advanced  in 
time. 


Efficioncy  (%) 


SIMULATION  OF  MIT  SHORT  CAVITY 


Beam  current  (A) 


A  :  Single  mode  simulations  with  fixed  a  (=1.93). 
+  :  Single  mode  simulations  with  varying  a. 

*  :  Multimode  simulations  with  varying  a. 

•  :  Experimental  measurements. 


Frequency  (GHz) 


EXPERIMENTAL  MEASUREMENT  OF 
SIDEBANDS  AT  MIT 


W.C.  Guss,  M.A.  Basten,  K.E.  Kreischer,  R.J.  Temkin. 
T.M.  AntGnsen  Jr.,  S.Y.  Cai,  G.  Saraph.  and  B.  Levush 


PREDICTED  AND  MEASURED  MODE  FREQUENCIES 


MARYLAND  LONG  CAVITY 


EFFICIENCY  VS.  CURRENT 
Maryland  Long  Cavity 


10  20  30  40  50 

Current 


Guss,  Basten,  Kreischer,  Temkin,  Antonsen.  Cai. 
Saraph,  and  Levush,  to  be  published  in  proceedings  of 
IRMW 


PARAMETRIC  MODE  COUPLING 
IN  HARMONIC  GYROTRONS 


Ph.  D.  Thesis 


Girish  Saraph 


PARAMETRIC  INSTABILITY 


Three  cavity  modes:  TEflpi,  TE(2P2,  and  TEUp2 

The  cyclotron  resonance  criteria  is  given  bv. 

^  -2  ^  U>3  Qo 

51  s2  >53  7 

The  resonance  criteria  for  the  parametric  instability  are 

as  follows: 

fa)  For  the  harmonic  numbers:  53  =  si  —  so. 

(b)  For  the  azimuthal  mode  indices:  /3  =  /1  j_ 

(c)  For  the  angular  frequencies:  u;3  ^ 

Two  characteristic  features  of  the  parametric  instability  are 
as  follows: 

(a)  Three  modes  are  phase-locked  together  , 

(b)  When  only  one  mode  has  high  amplitude  then  the 
small  signal  gains  of  the  other  two  modes  match. 

When  the  three  modes  .are  parametrically  coupled,  then  differ¬ 
ent  types  of  mode  interactions  can  take  place,  viz.  parametric 
excitation,  decay  of  one  or  two  modes,  co-existence  of  three 
modes,  or  cyclic  mode  hopping  between  the  three  modes. 

Notation  used: 

Moae  1  .  TE(lPl  at  the  fundamental  frequency  =  1) 

Mode  2  :  TEi2P2  at  the  second  harmonic  frequency  (s2  =2) 
Mode.  3  :  TEi3Pz  at  the  third  harmonic  frequency  (s3  =3) 


Normalized  Mode  Amplitude 


Uilll 


Normallz 


6 

alized  Time 


REGION  OF  PARAMETRIC  EXCITATION  OF 


A  THIRD  HARMONIC  AT  94  GHZ 


2.0  I""' 11  1,11 
"  M  =  7-5 


Unstable  fieaion 


DESIGN  OF  A  THIRD  HARMONIC  GYROTRON 

AT  94  GHZ 


Moae  I  :  TE(Z.l)  at  si  =  i 
Mode  2  :  TE(  1.3)  at  s2  =  2 
Mode  3  :  177(4.3)  at  s3  =  3 

Selected  Quantities: 

Wall  Radius  rw  =0.65  cm. 

Beam  Radius  rh  =  0.42  cm. 

Quality  Factors:  Q1  =  300. 

Q2  =  1200.  Q3  =  2700 

Normalized  Parameters: 

(j.  =  7.5.  tv  =  86.  Il  =  230. 

<5  =  0-7*  e.3  —  3.18.  r/j_  =  0.412 

Corresponding  Physical  Quantities: 

Interaction  Length  Lc  =  4.14  cm 
Operating  Frequency  fz  =  94  GHz 
Magnetic  Field  £0  =  1.15  T 
Beam  Voltage  V  =  31  A  V 
Beam  Current  I  =  15.75  ,4mr> 

A 

Pitch  Angie  a  =  1.9 
Net  Efficiency,  rj  =  32.3  % 

Total  Power  Output  P()  =  15 8  kW 


Multi-frequency  simulation  of  high 
gain  fel  oscillators 


Radiation 


+  vg^]  5a(z,t)  = 


7rl  VgVw(z)  ev  'er*  C 

U  k0  Aeff  vz 


<exp(-i  \| /)> 


Electrons 


3  3  Icq 

+  Vz3z^  ^ 1  T  vw(z)  ey  -er  C  exp(  i  \j/)  +  c.c. 


3  3 

f3t  +  Vz3?  v  =  (k  +  kw)vz-co 


Proposed  FEM  design 


Beam  voltage 

Beam  current 

Beam  emittance 

Beam  radius 

Wiggler  period 

Peak  field  wiggler  section  1 

Peak  field  wiggler  section  2 

No.  periods  wiggler  section  1 

No.  periods  wiggler  section  2 

Inter  wiggler  gap 

Waveguide  mode 

Waveguide  width 

Waveguide  height 

Cutoff  frequency 

Cavity  length 

Power  reflection  coefficient 


1.75  MV 
12  A 
0 

.15  cm 
4  cm 
2kG 
1.7  kG 

25  (23  full  periods) 
19  (17  full  periods) 
6  cm 

HE ii  rectangular 

1.5  cm 

2.0  cm 

12.49  GHz 

382  cm 

.21 


GAIN 


Output  Spectrum  versus  Time 


frequency 


Modeling  of  harmonic  gyroklystrons, 
phase  locked  gyro-oscillators, 
and  gyro-twystrons 


Recent  Topics 

1 .  Theory  of  the  relativistic  gyrotwystron 

input:  cavity 

output:  travelling  wave  amplifier 

2 .  Two  harmonic  prebunching  of  electrons  in 
multi-cavity  gyrodevices 

input  cavity:  fundamental 
second  cavity:  fundamental/second 
harmonic 

output  cavity:  fundamental- fourth  harmonic 

3 .  Theory  of  phase-locked  gyrotrons 
operating  at  cyclotron  harmonics 

input  cavity:  fundamental 
output  cavity:  second  harmonic 


3.  Theory  of  phase-locked  gyrotrons 
operating  at  cyclotron  harmonics 

Locking  Bandwidth  versus  Beam  Current 


cold  beam 


Frequency  Shift 


3.  Theory  of  phase-locked  gyrotrons 
operating  at  cyclotron  harmonics 


Locking  Bandwidth  versus  Beam  Current 


5  a/a  =  .2 


Frequency  Shift 


3.  Theory  of  phase-locked  gyrotrons 
operating  at  cyclotron  harmonics 

Locking  Bandwidth  versus  Beam  Current 

5a/a  =  .2  with  compensation 


Frequency  Shift 


Current  Efforts 

Adapt  individual  models  into  a  flexible  simulation 
code  for  design  and  analysis  of 

Phase-locked  harmonic  gyrotron  oscillators 

Gyro-klystrons 

Gyro-twystrons 

Phigtrons 

Included  effects 

Frequency  multiplication 

Multiple  sections 

Backward  waves 


NUMERICAL  SIMULATION 
OF 

SLOW  WAVE  DEVICES 


MODELING  OF 

vacuum  backward  wave  oscillators 
plasma  filled  backward  wave  oscillators 


EXAMPLES 

homogeneous  slow  wave  structure 
effect  of  the  finite  pulse  duration 
cyclotron  absorption  effect 
effect  of  the  non-synchronous  harmonic 
effect  of  plasma 


non-homogeneous  slow  wave  structure 
efficiency  enhancement 


THEORETICAL  MODEL 


*  representation  of  the  field 


^  ,  i(k  z-cot )  i(k  z-cot ) 

E(x,t)={s  (■ z,t)e  ~  E  (x,k  )+£(z,t)e  +  E  (x,k  )}+c.c 
~  P  ~  +  P  + 


i(k  z-cot )  i(k  z-cot ) 

B{x,i)={£  (. z,t)e  B  (x,k  )+£(z,t)e  +  B  (x,k  )}+c.c 

-  p  -  +  p  + 


£ 

+ 


(z,0— >slowly  varying  envelope 


7 


E  ( x,k  ) 

p  ± 

periodic  functions 

B  (x9k  ) 

P  ± 


ink  z  A  ^  2ir 

E  (x,k)=lE  (x.k  )e  0  *■>'  J 

p  ±  n  1  ± 


S  {x,k  )=S  5  (x  ,k  )e‘nk°Z 
p  ±  n  ±  ± 


g 


de_  de_  rd  . 

~df + v°'-^  =  ~^dz 


f-T,u  dt  I  <i2x±E;(x,  k_) .  j e-(*— "')/[/,  (9) 


electromagnetic  energy  per  Ie2|  contained  in  one 
period 


U  =  Jad  dz  J  d2xJ^L±^p\2 


(10) 


group  velocity  for  the  vacuum  backward 
electromagnetic  wave 


VS-  =  Jo  dzJd2x±^-zo  ■  (Ep(x,fc_) 


x  B!(x,  k_))/u 


(11) 


5 


.  particle  motion 


dp 

dt 


q 


E  +  lxlB  +B 


ext 


A 


dt 

•  assume  that  only  two  spatial  harmonics  of  the 
backward  wave  can  interact  strongly  with  the 
beam 


i(k  z—iot) 
E  (x,t)  =  8  e  ~ 


-i{k  z) 

ao  Eo(x,k-  +  a-iE-i^x,k-  6  ° 


+  c.c 


B  [x,t)  =  £  C  ~ 


i(k  z-coO 


a  B  (x,k 
0  0  - 


)  +  a  (. x,kj  e 


-i(k  z) 
0 


+  c.c 


INTERACTION  EFFICIENCY  [%]  INTERACTION  FREQUENCY  (GHz) 
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infinitely  long 
structure  calculation 


o  o  o 
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o  theory 
*  experiment 


Double  Cathode 
AKJ/0=  48/51.5  mm 
B=  14.8  kG 


300  400  50 

ELECTION  GUN  VOLTAGE  (kV) 


*  experiment 


300  400  500  600 

ELECTION  GUN  VOLTAGE  (kV) 


Cyclotron  Absorption  Effect 


B(kG) 


I 

( 


Only  th 


O 


t<nsec> 


synchronous  harmonic 


Modeling  of  plasma  filled  Backward  Wave  Oscillator 

Ph.  D.  Thesis 
Susan  Miller 


Efficiency  enhancement 
of  relativisitc  backward 
wave  oscillator 

©  vary  coupling  impedence 
along  the  interaction  region 


9  sudden  increase  in  phase 
velocity  of  the  synchronous 
wave. 


Efficiency  enhancement 
of  high  power  relativistic 
backward  wave  oscillator 


Homo ^ e  ;S”frn 


’Vc**-  pi*) 


Efficiency  vs.  time  (500  keV  &  100  A  beam) 
with  space  charge  neglected  for  uniform 


nromalized  |Ez| 


Amplitude  of  synchronus  field  (red)  and 
space  charge  harmonics  fields  (black) 
vs-  axial  distance 

(500  keV  &  5  kA  uniform  sturcture) 


ay 

5i  = 


"radiation  term" 


-iIWRe{|1R„ein¥<e'inv>}’ 


particle  phase:  \j /  =  k_z  -  cot 


Rn  =  geometric  space  charge  reduction  factor 


MAGIC  simulation  of  S  —  band  BWO 
with  uniform  structure 
(12  and  14  periods  long) 
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MAGIC  simulation  of  S-band  I3W0 
with  non-uniform  structure 
(transition  at  3.5  and  <1.0  periods) 
Output  power  vs.  time 
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ns 


Current  efforts 


PIC  modeling  of  non-homogenous  structures 
modeling  of  plasma  filled  backward  wave  oscillator 


Possible  new  efforts 


adopting  our  models  to  simulate  multi-frequency  phenomenon 
in  slow  wave  devices  such  as  TWT’s 

1 .  self-excitation  of  spurious  oscillations 

2.  simultaneous  amplification  of  a  number 
of  signals  at  different  frequencies 


